Universal trend of the information entropy of a fermion in a mean field 
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We calculate the information entropy of single-particle states in position-space S r and momentum- 
space Sk for a nucleon in a nucleus, a A particle in a hypernucleus and an electron in an atomic 
cluster. It is seen that S r and Sk obey the same approximate functional form as functions of the 
number of particles, S r (or Sk) = a + bN 1 ^ 3 in all of the above many-body systems in position- 
and momentum- space separately. The net information content S r + Sk is a slowly varying function 
of iV of the same form as above. The entropy sum S r + Sk is invariant to uniform scaling of 
coordinates and a characteristic of the single-particle states of a specific system. The order of 
single-particle states according to S r + Sk is the same as their classification according to energy 
keeping the quantum number n constant. The spin-orbit splitting is reproduced correctly. It is also 
seen that S r + Sk enhances with excitation of a fermion in a quantum-mechanical system. Finally, 
we establish a relationship of S r + Sk with the energy of the corresponding single-particle state i.e. 
S r + Sk = fcln(/i_E + v). This relation holds for all the systems under consideration. 
PACS: 89.70.+c; 36.40.+d; 31.10.+z; 21.60.-n 
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I. INTRODUCTION 

The information entropy for a continuous probability distribution p(x) in one dimension is defined by the expression: 



^: s = - 1 p{x)H P {x))dx (i) 

. where Jp(x)dx = 1. 

S is measured in bits if the base of the logarithm is 2 and nats (natural units of information) if the logarithm is 
natural. It represents the information content of a probability distribution as well as a measure of uncertainty of the 
1 corresponding state. We note that the information and thermodynamic entropy are different concepts but can be 
| connected employing some assumptions. 

An important step in the past was the discovery in [1] of an entropic uncertainty relation (EUR), which for a 
i— H ' three-dimensional system has the form: 

o : 

S r + S k >3{l + In tt) ^ 6.434, (h = 1) (2) 

where 

,s', = -/,.(,•) In ,.(,-Wr (>!) 

is the information entropy in position-space 

S k = - I n(k)lnn(k)dk (4) 



is the information entropy in momentum-space and p(r), n(k) are the density distributions in position- and momentum- 
space respectively, normalized to unity. 

The lower bound in (Q) is attained for gaussian density distributions. The physical meaning of the above inequality 
is the following: an increase of Sk corresponds to a decrease of S r and vice versa, which indicates that a diffuse density 
distribution n(k) in momentum space is associated with a localised density distribution p(r) in configuration space 
and vice versa. 

Relation (2) represents a strengthened version of Heisenberg's uncertainty principle for two reasons: first EUR leads 
to Heisenberg's uncertainty relation but the inverse is not true. Second the right-hand-side of EUR does not depend 
on the state of the system, while in Heisenberg's relation does depend. It is obvious from (^) and (|J) that S r , Sk 
depend on the unit of length in measuring p(f) and n(k). However, the important quantity is the entropy sum S r + Sk 
(net information content of the state) which is invariant to uniform scaling of coordinates. 
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Information entropy was employed in the past for the study of quantum mechanical systems [1-9]. Recently jlO| we 
studied the position- and momentum- space information entropies S r , Sk respectively for various systems: the nuclear 
density distribution of nuclei, the electron density distribution of atoms and the valence electron density distribution 
of atomic clusters. We showed that a similar functional form S — a + b In TV for the entropy as function of the number 
of particles N holds approximately for the above systems. We conjectured that this is a universal property of a 
many-fcrmion system in a mean field. 

The concept of information entropy proved also to be fruitful in a different context jj| . We used the formalism of 
Ghosh, Berkowitz and Parr |pd[| within the ground state density functional framework, to define the concept of an 
information entropy associated with the density distribution of a nuclear system. It turned out that S increases with 
the quality of the wave function and can serve as a criterion of the quality of a nuclear model. 

Another interesting result |L2j is the fact that the entropy of an A-photon state subjected to Gaussian noise increases 
linearly with the logarithm of A. 

Encouraged by previous work we attempt in the present paper to calculate S r , Sk for the wave functions of single- 
particle states (instead of the total densities as in |n|) for various systems i.e. a nucleon in a nucleus, an electron in 
an atomic cluster and a A particle in a hypernucleus. We employ for these systems models existing in the literature. 

Our aim is to investigate the dependence of S r , Sk on the excitation of a fermion in a quantum- mechanical system 
as well as its dependence on the system under consideration and the number of the particles N . We also attempt to 
connect the information entropy with the energy of the single-particle state. The study of the dependence of S on the 
quantum state of a system is also interesting (as stated in M) for two reasons: (i) The information-theoretical and 
physical entropy are connected via Boltzmann's constant kg by the Jayne's relation S p h ys — ksSmf. Thus one can 
ascribe to any quantum object a certain value of its physical entropy S p h ys if one calculates Si n f. (ii) It is interesting 
to know the value of the information entropy which is a measure of the spatial " spreading out" of the wave function 
for various states of various systems. 

The present paper is organized as follows: In sec. 2 we calculate S r , Sk for single-particle states of a nucleon in 
a nucleus as function of the number of nucleons A using the simplest model available i.e. the harmonic oscillator 
potential and a more realistic one (Skyrme). In sec. 3 we calculate S r , Sk for a A particle in a hypernucleus employing 
a simple and (semi-) analytical relativistic model. In sec 4 we determine S r , Sk for the single-particle states of an 
electron in atomic (metallic) clusters using the Woods-Saxon potential. In sec. 5 we present a relationship of S r + Sk 
with the energy. Finally, sec. 6 contains a discussion of our results (comparison of sec. 2, 3 and 4) and our conclusions. 



II. INFORMATION ENTROPY FOR A NUCLEON IN A NUCLEUS 

The information entropy S r in position-space for a single-particle wave function ijj(r) is defined as 

S r = -J\ij(r)\ 2 \n\^(r)\ 2 dr (5) 
while the entropy Sk in momentum-space is 

S k = - /|0(fc)| 2 ln|0(fc)| 2 dfc (6) 



where <p(k) is the Fourier transform of tp{r). 

In this section we calculate S r and Sk for the single-particle states Is, lp, Id, • • • of a nucleon in a nucleus in 
the framework of the harmonic oscillator (HO) model. We use for the HO parameter the well-known expression 
huu = 41A- 1 / 3 MeV. 

We find that the value of hut is important only for S r , Sk, while the net information content S — S r + Sk is 
independent of hui and consequently of A. It depends only on the state under consideration and characterises it. 
These values for the states Is, lp, Id and 2s are 6.4341, 7.8388, 8.6651 and 8.3015 respectively. 

However, the HO model is a simplification. Thus, we employed a more realistic parametrization of the nuclear 



mean field i.e. the Skyrme (Sk III) interaction 
We choose to work with protons. However, simi 



[13J. In this model protons and neutrons move in different potentials, 
ar results can be obtained for neutrons. We found that the values for 
S r , Sk obtained from the wave functions of single-particle states calculated according to Sk III are represented well 
by the expression 

S r (orS k )=a + bN 1 / 3 (7) 

while S r + Sk is a slowly varying function of A of the same form as (Q). The values of the parameters are shown in 
Table 1. 
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In Fig. la we plot our fitted expressions (Sk III) S r (or Sk) = a + bN 1 / 3 for the entropies S r , Sk, S r + Sk of ls-states 
as functions of N 1 ^ 3 . The lines correspond to our fitted expressions, while the corresponding values of our numerical 
calculations are denoted by squares for S r , circles for Sk and triangles for S r + Sk- Similar graphs can be plotted for 
the higher states lp, Id, 2s, ■ ■ ■. From Fig. la we see that the values of the entropies are represented well by our fitted 
expressions. In Fig. lb we compare the sum S r + Sk = a + WV 1 / 3 for various single-particle states. We observe that 
the entropy sum S r + Sk enhances with the excitation of the single-particle states. We see that S r + Sk is a slowly 
varying function of N. We also note that the spin-orbit splitting is reproduced correctly i.e. the state lpa/2 is lower 
than lpi/2 e.t.c. (as for the energy) although their difference is small and cannot be shown in the figure. 

III. INFORMATION ENTROPY FOR A A IN A HYPERNUCLEUS 



We employ a simple and (semi-) analytical relativistic model of a hypernucleus of ref. jbljjlfj], where the Dirac 
equation with a scalar potential Us(r) and the fourth component of a vector potential Uy{f) was considered in the 
case of rectangular shapes of these potentials with the same radius: 

R = r Ali 3 re 

In |Q the Dirac equation was solved and gave the wave functions G(r) and F(r) for the large and small components 
for a A particle in a hypernucleus. These components can be found in relations @ and (0) of Q. 
The Dirac spinors in terms of large (G) and small (F) components can be expressed: 

iG{r)/r 



* - I F(r)/r ) ^ 



The density distribution of a A in position space is: 
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and the normalization is 



P('-) = S [G 2 (r)A 2 +F 2 (r)/r 2 ] (9) 



47r / p(r)r dr = 1 
J o 

In momentum-space we have: 

"' M ~ . Y{k) 

where X(k) and Y(k) are the Fourier transforms of G(r)/r and F(r)/r respectively. Thus the density distribution in 
momentum-space is given by: 

n(fc) = i-[X 2 (fc) + y 2 (fc)] (11) 

and the normalization is: 

/CO 
n{k)k 2 dk = 1 
o 

The information entropies of the A particle are calculated according to the relations: 

S r = -Air J p{r) In p{r)r 2 dr (12) 

S k = -4tt / n(k) In n(k)k 2 dk (13) 



where p(r) and n(k) are given by (O) and (11) respectively 
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For the depths of the potential we used the values (h): D+ = 30.55 MeV, L>_ = 300 MeV, r = 1.01 fm and 

1 /3 

the radius parameter R = r^Acore obtained by fitting the experimental binding energies of the ground state of the A 
particle. In the following we put A core — N — number of particles. 

Next we fitted the expressions S r (or Sk) = a + bN 1 / 3 to the values of S r , Sk calculated from ( |l2| ) and (|l^) and 
found that these values are represented well. The values of the parameters a and b for various states are shown in 
Table 1. 

In Fig. 2a we plot our fitted expressions for S r , Sk, S r + Sk as functions of N 1 / 3 for the Is state. This is done for 
a A in a hypernucleus in a similar way as for a nucleon in a nucleus (Fig la). Similar graphs can be plotted for the 
higher states. In fig. 2b we compare the sum S r + Sk for various single-particle states of a A (similar with Fig lb). 
The spin orbit splitting is reproduced correctly as in nuclei (Sec. 2). 



IV. INFORMATION ENTROPY FOR AN ELECTRON IN AN ATOMIC CLUSTER. 



We consider atomic (metallic) clusters composed of neutral sodium atoms, where the electrons move in an effective 
radial electronic potential parametrized by Woods-Saxon potential of the form: 

—Vn 

1 + cxp[(r — R)/a\ 

with Vq = 6 eV, R = tqN 1 / 3 , Tq = 2.25 A and a = 0.74 A. For a detailed study regarding the parametrization of 
Ekardt's potentials see ref |T^| . 

We found the wave functions of the single-particle states in configuration space and by Fourier transform the 
corresponding ones in momentum space by solving numerically the Schrondinger equation for atomic clusters for 
various values of the number of valence electrons N . Using the above wave functions, we calculated the information 
entropies S r , Sk (relations (0) and (0) ) for the single-particle states instead of the total density distributions as in 
ref. |lC[ |. Then we fitted the form S r (or Sk) — a + bN 1 / 3 to these values and found that these expressions represent 
well the values of S r ,Sk- In Fig. 3a we plot S r , Sk and S r + Sk as functions of N 1 / 3 (similar as Fig. la, 2a) and in 
Fig. 3b we compare S r + Sk for various states (similar as in Fig. lb, 2b). In Table 1 we present the values of the 
parameters a and b which were obtained from the fitting. 



V. RELATIONSHIP OF THE INFORMATION ENTROPY WITH THE ENERGY OF SINGLE-PARTICLE 

STATES 

In Fig. 4 we plot S r + Sk, obtained with the HO model of the nucleus, versus the energy of the single-particle 
states. We use hu; = 41 A~ x / 3 with A — 208 (Pb) and keep the quantum number n equal to 1. A fitting procedure 
gives for n = 1, the relation: 

S = k\n(jiE + u) (15) 

where k = 2.0206, fi = 3.5373 MeV -1 and v = -12.5320. Similar relations hold for n > 1. 

Next we plot the sum S r + Sk as function of the energy E of single-particle states for a proton in a nucleus according 
to Sk III interaction for 208 Pb (Fig. 5) and an electron in atomic cluster with N = 198 (Fig. 6) for n = 1. Similar 
curves hold for higher values of n > 1. In both cases the dependence of S r + Sk on E can be represented well by the 
functional form (jl5|). The values of the constants are the following 

k = 1.5262, n = 17.3043 MeV -1 , = 793.109 for a proton in a nucleus 



k = 1.2386, fj, = 1481.48 eV~\ f = 8730.52 for an electron in a cluster 

A similar relation may be obtained for a A in a hypernucleus but the number of values of S r + Sk available is 
small. It is the first time in our research on information entropy that we observe such a relationship of S r + Sk with 
a fundamental quantity as the energy. We note that in Fig. 5 there are pairs of points with almost the same S r + Sk 
which reproduce the spin-orbit splitting. 
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VI. DISCUSSION AND CONCLUSIONS 



Comparing our results in sections 2, 3, and 4, we see that a similar functional form S r (or Sk) = a + bN 1 / 3 describes 
well the information entropies S r , Sk of the single-particle states for a nucleon in a nucleus, a A in a hypernucleus 
and a valence electron in an atomic cluster, although the single-particle potentials are different. We conjecture that 
this is a universal trend of the information entropies S r , Sk for a fermion in a mean field, while the net information 
content S r + Sk of the single-particle states of a fermion in a mean field is a slowly varying function of N of the form 
S — a + bN 1 / 3 for the systems considered above. For nuclei and the simple HO potential S r + Sk is exactly a constant 
independent of N i.e. 6 = 0. 

We note that in E0[ we found the universal property S = a + b In N for the total density distributions of various 
systems. 

In both cases it is not clear why S depends linearly on IniV (total densities) or linearly on TV 1 / 3 (single-particle 
states) but we note that in atomic physics there is already a connection of the information entropy with experiment 
i.e. with fundamental and/or experimental quantities e.g. the kinetic energy or the magnetic susceptibility. Both 
characteristics have been used in the study of the dynamics of atomic and molecular systems |Tij] . This connection 
established the information entropy as an interesting entity for atomic physics. In the present paper we obtained a 
relationship of S r + Sk with a fundamental quantity as the energy of the single-particle states, i.e. S = k ln(/^_B + v). 
It is remarkable that the same functional form holds for various systems. 
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FIG. 1. (a) Values of the information entropies S r (squares), Sk (circles) and S r + Sk (triangles), calculated numerically, 
versus the number of particles N. These values correspond to the single-particle states of a proton in various nuclei, according 
to the Sk III interaction. The lines correspond to our fitted expressions S r (or Sk) = a + bN 1 ^ 3 . (b) Comparison of the sum 
S r + Sk for various proton single-particle states. 




FIG. 2. The same as in Fig. 1 for a A in hypernuclei employing a relativistic model. 
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FIG. 3. The same as in Fig. 1 for an electron in atomic clusters with the Woods-Saxon potential. 
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FIG. 4. The values of the energy E (squares) of single-particle states for a nucleon in Pb according to the HO model for 
n = l. The line corresponds to our fitted expression S = k \n(nE + v). 
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FIG. 5. The same as in Fig. 4 but for a proton in Pb 208 according to the Sk III interaction for n = 1. 
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FIG. 6. The same as in Fig. 4 but for an electron in atomic clusters (N=198) using the Woods-Saxon potential for n = 1. 
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TABLE I. Values of the parameters a and b which appear in the expressions S r (or Sk) = a + WV 1 / 3 for a nucleon (proton) 
in nuclei according to the Sk III interaction, a A in hypernuclei according to a relativistic model and an electron in atomic 
clusters with Woods-Saxon potential. 



Case 


State 


Sr 

a 


b 


a 


s k 


b 


Sr ~t~ Sk 

a 


b 


Nucleus 


lSl/2 


3.0831 


0.8652 


3.2353 




-0.8140 


6.3217 


0.0501 






4.2824 


0.6368 


3.6256 




-0.6688 


7.9084 


-0.0322 




lPl/2 


4.2724 


0.6235 


3.6799 




-0.6675 


7.9521 


-0.0439 




1^5/2 


4.7500 


0.5743 


4.0980 




-0.6513 


8.8480 


-0.0771 




1^3/2 


4.9553 


0.5042 


4.0080 




-0.6071 


8.9618 


-0.1024 




2*1/2 


5.2456 


0.3364 


3.4355 




-0.4641 


8.6756 


-0.1260 


Hyper- 


lSl/2 


3.5817 


0.4967 


2.8303 




-0.4756 


6.4120 


0.0214 


nucleus 


1P3/2 


4.4347 


0.3789 


3.5123 




-0.4021 


7.9475 


-0.0232 




± PL/2 


4.3764 


3835 


3.6575 




-0.4199 


8.0342 


-0.0364 




1^5/2 


5.2553 


0.2462 


3.6249 




-0.2938 


8.8803 


-0.0475 




1^3/2 


4.9910 


0.2819 


4.0503 




-0.3543 


9.0414 


-0.0724 


Cluster 


Is 


4.3038 


0.7113 


2.0923 




-0.6883 


6.3960 


0.0232 




lp 


5.1114 


0.6135 


2.7700 




-0.6299 


7.8816 


-0.0163 




Id 


5.5191 


0.5636 


3.1420 




-0.5842 


8.6611 


-0.0205 




2s 


5.3858 


0.4918 


2.9672 




-0.5301 


8.3536 


-0.0383 
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